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Abstract. The major problem in the cosmological nucleosynthesis is the evaluation of the re- 
I— — I action rate. The present scenario is that the standard thermonuclear function in the Maxwell- 

Boltzmann form is evaluated by using various techniques.The Maxwell-Boltzmannian ap- 
dj proach to nuclear reaction rate theory is extended to cover Tsallis statistics (Tsallis, 1988) 

^ and more general cases of distribution functions. The main purpose of this paper is to in- 

vestigate in some more detail the extended reaction probability integral in the equilibrium 
thermodynamic argument and in the cut-off case. The extended reaction probability integrals 
^ will be evaluated in closed form for all convenient values of the parameter by means of residue 

^ calculus. A comparison of the standard reaction probability integrals with the extended re- 

^ action probability integrals is also done. 
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lo 1 Introduction 

o 

The evolution of universe is due to the thermonuclear reactions which are taking 
^ place in hot cosmic plasma. The main concept behind the description of cosmic nucle- 

Q osynthesis is the rate of nuclear reactions synthesizing light nuclei into heavier ones. If 

^ we closely examine the nuclear cross section theoretically and experimentally, we can 

find the analytical representations of thermonuclear reactions. Many researchers were 
. ^ looking for these analytic representations of the reaction rate probability integrals for 

^ the last few decades. Many approaches have been made on the study of thermonu- 

d clear reactions(Haubold and John, 1978; Haubold and Mathai, 1985, 1998; Anderson 

et al, 1994, Saxena et al, 2004). These studies will be effective only when the reaction 
probability integrals are expressed in computable series representations(Mathai and 
Haubold, 1988; Haubold and John, 1982). The derivations of closed-form representa- 
tions of nuclear reaction rates and the approximations on them are based on the theory 
of generalized special functions. 
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In the production of neutrinos in the gravitationally stabihzed solar fusion reactor, 
due to the memory effects and long-range forces a possible deviation of the velocity 
distribution of plasma particles from Maxwell-Boltzmann is noted(Coraddu et al, 1999; 
Lavagno and Quarati, 2002; Lavagno and Quarati , 2006). This was initiated by Tsallis' 
non-additive generalization of Boltzmann-Gibbs statistical Mechanics. Tsallis statistics 
covers Boltzmann-Gibbs statistics( Tsallis, 1988; Gell-Mann and Tsallis, 2004; Tsallis, 
2004). The extension of the nuclear reaction rate theory from Maxwell-Boltzmann the- 
ory to Tsallis theory was done by Saxena et al (2004), Mathai (2005), Mathai and 
Haubold (2007). In this scenario Mathai introduced a more general distribution func- 
tion which can be incorporated in the reaction rate theory by appealing to cntropic 
and distributional pathways. If Mathai's pathway model (Mathai, 2005; Mathai and 
Haubold, 2007) is introduced in the reaction probability integrals one goes into a wider 
class of integrals where the standard reaction probability integrals becomes a limiting 
case. 

The paper is organized in the following way. Section 2 contains the standard rep- 
resentations of the non-resonant thermonuclear reaction rates. In section 3 we give 
the basic definitions that we use in our discussion. Section 4 gives an outline of the 
extensions of the reaction probability integrals using the pathway model and also es- 
tablish the series representations of the extended integrals Iia and ■ Section 5 gives 
a comparison of the extended integrals with the standard integrals. 

2 Standard representations of non-resonant 
thermonuclear reaction rates 



For the evaluation of the reaction rate of the interacting particles i and j we 
have to consider the energies distributed between the particles. The reacting particles 
in the astrophysical plasma follows a Maxwell-Boltzmann distribution. From Mathai 
and Haubold (1988) we can see that the expression for the reaction rate rij of the 
reacting particles in the non-degenerate environment is 

Ea{E)e-^dE (2.1) 

= niUj {(71') 

where rii and rij are the number densities of the reacting particles i and j, the reduced 
mass of the particles is denoted by /j, — , T is the temperature, k is the Boltzmann 

constant, the reaction cross section is cr{E) and the kinetic energy of the particles in 

2 

the center of mass system is £■ = ^ where u is the relative velocity of the interacting 
particles i and j. 

We write (au) to indicate that it is an appropriate average of the product of the 
cross section and relative velocity of the interacting particles. For detailed physical 
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reasons see Haubold and Mathai (1984, 1986). 



2.1 Standard non- resonant thermonuclear function 

When two nuclei of charges Zi and Zj are coUiding at low energies below the coulomb 
barrier, the reaction cross section for the non-resonant nuclear reactions have the form 
(Haubold and Mathai, 1998; Bergstroem et al, 1999, Mathai and Haubold, 2002) 

a^E) = ^^e-^^''^^) (2.2) 



with 



E 

V{E) = (f)^^ (2.3) 



where r){E) is the Sommerfeld parameter, h is the Planck's quantum of action, e is 

the quantum of electric charge, the cross section factor S{E) is often found to be a 
constant or a slowly varying function of energy over a limited range of energy (Mathai 
and Haubold, 1988). The cross section factor S{E) can be parameterized by expanding 
in terms of the power series about the zero energy because of its slow energy dependence. 
S{E) can be expressed as 

S{E) - S{0) + —E + i^-^E , (2.4) 

where S(0) is the value of S{E) at zero energy and 5"(0) and ^"'(O) are the first and 
second order derivatives of S{E) with respect to energy evaluated at £■ = 0, respectively. 
Then 



i/=0 ^ ^ 

1 2 

where x — and b = (2^) ^ ^'^^^ . The standard case of the thermonuclear function 
contains the nuclear cross section cr(£'), the energy dependent cross section factor S{E) 
and the steady-state Maxwell-Boltzmann distribution function. 

The collision probability integral for non-resonant thermonuclear reactions in the 
Maxwell-Boltzmannian form is (Haubold and Mathai, 1984) 

hiu,l,b,-) = x'^e— ^-"Mx. (2.6) 

We will consider the general integral 

poo 

/i(7- l,a,6,p) = / x^-V'^^-^^^'dx, a>0,6>0,p>0. (2.7) 
Jo 
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2.2 Non-resonant thermonuclear function with high energy 
cut-off 



It is assumed that the thermodynamic fusion plasma is in exact thermodynamic 
equihbrium. But the cut-off of the high energy tail of the Maxwell-Boltzmann distri- 
bution function in (2.6) results in a modification of the closed form representation of 
the appropriate quantity {av) which is given by 

1, h,\)= ! x'^e-^-^^'Mx, 6 > 0, d<oo. (2.8) 
Again we consider the general form of the integral (2.8) as 



ir (7 - 1, a, 6, p) = / a;^-^e-"^-''^"''da;, a > 0, 6 > 0, p > 0, < oo. (2.9) 



Jo 

For physical reasons for the cut-off modification of the Maxwell-Boltzmann distri- 
bution function of the relative kinetic energy of the reacting particles refer to the paper 
Haubold and Haubold and Mathai (1984). 

2.3 Modified non-resonant thermonuclear function with de- 
pleted tail 

A depletion of the high energy tail of the Maxwell-Boltzmann distribution func- 
tion of the relative kinetic energies of the nuclei in the fusion plasma is explained in 
Haubold and Mathai(1986); Haubold and John (1982); Kaniadakis et al (1997,1998). 
The ad hoc modification of the Maxwell-Boltzmann distribution for the evaluation of 
the non-resonant thermonuclear reaction looks like a depletion of the high energy tail 
of the Maxwell-Boltzmann distribution. If their exists a possibility of such a modifica- 
tion a remarkable change of the views of astrophysical nucleosynthesis and controlled 
thermonuclear fusion may arise. 

The integral form [av) in comparison with strict Maxwell- Boltzmannian case, we 
have the integral 




6 > 0. 



(2.10) 



We will consider the general integral of the type 




oo 



X 



,'y—l —ax —bx ^ 



dx. 



(2.11) 



where 2; > 0, a > 0, 6 > 0, p > 0. 



3 Mathematical preliminaries 
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The basic quantities which we need in our discussion will be given here. The gamma 
function denoted by T{z) for complex number z is defined as 



r(^) = r 

Jo 



t'-^e-*dt,^{z) > (3.1) 



where 3?(-) denotes the real part of (•). In general r{z) exists for all values of z, positive 

or negative, except at the points z = 0, —1, —2, ■ ■ ■ . These are the poles of T{z). But 
the integral representation holds for the real part of z to be positive. Another important 
result that wc use in our discussion is the psi function. The psi function which is denoted 
by tp{z) is the logarithmic derivative of a gamma function and is defined as 

J d[r(4] 
V'(^) = ^[lnr(^)] = or Inr(z) = i^{x)dx. (3.2) 

One property of the psi function that we will use is 

^|;{l + n) = l + l + \ + ■■■ + --^ (3.3) 
2 6 n 

where 7 is the Euler's constant and 7 = 0.5772156649 • • • , n = 1, 2, 3, • • • . The G- 
function which is originally due to C. S. Meijer in 1936 (See Mathai, 1993; Mathai and 
Saxena, 1973) is defined as a Mellin-Barnes type integral as follows: 



where i = y — 1, L is a suitable contour and z 0, m, n,p, q are integers, < m < q and 
< n <p, the empty product is interpreted as unity and the parameters oi, 02, • • • , Oj, 
and bi, &2, • " " > ^-re complex numbers such that no poles of r{bj + s), j — 1, • ■ • ,m 
coincides with any pole of r(l — Uk — s), k — 1, • • ■ ,n; 

—bj — u ^ 1 — Qk + X, j — 1, ■ ■ ■ ,m; k — 1, ■ ■ ■ ,n; u, X — 0,1, ■ ■ ■ . 

This means that — hj 7^ z/ + A + 1 or a^. — hj is not a positive integer for j = 
1, • • • , m; k = 1, • • • ,n. We also require that there is a strip in the complex s-plane 
which separates the poles of T{bj + s), j — 1, - • ■ ,m from those of r(l — — s), k — 
1, • • • ,n (For the existance conditions and properties of G-functions see Mathai (1993)). 

Next we need the pathway model of Mathai (2005). When fitting a model to ex- 
perimental data very often one needs a model with a thicker or thinner tail than the 
ones available from a given parametric family, or sometimes we may have a situation 
of the right tail cut-off. In order to take care of these situations and going from one 
functional form to another, a pathway parameter is introduced, see Mathai (2005) and 
Mathai and Haubold (2007). By this model we can proceed from a generalized type-1 
beta model to a generalized type-2 beta model to a generalized gamma model when the 



5 



variable is restricted to be positive. For the real scalar case the pathway model is the 
following: 

f{x) ^c\x\'^-^[l-a{l-a)\xf]^, a>0,S> 0, l-a(l-Q;)|x|'^ > 0,7 > 0,7^ > (3.5) 

where c is the normalizing constant and a is the pathway parameter. When a < 1 the 
model becomes a generahzed type-1 beta model in the real case. This is a model with 
the right tail cut-off. When a > 1 we have 1 — « = — (a — 1),q;>1so that 

f{x) = c|xp-^[l + a{a - 1)1x1'']-^, (3.6) 

which is a generahzed type-2 beta model for real x. When a — > 1 the above 2 forms 
will reduce to 

f{x) ^clxp-'e-"'"''. (3.7) 
Observe that the normalizing constant c appearing in (3.5), (3.6) and (3.7) are different. 

4 Extended thermonuclear function through 

pathway model 

When Mathai's pathway model is introduced in the reaction probability integrals 
we get a wider class of integrals. Then the standard reaction probability integrals 
become particular cases of the new family of integrals. Through the pathway parameter 
a we move to a wider class of integrals as a — )■ 1 we get the reaction rate probability 
integrals. If Maxwell-Boltzmann is the stable situation , many unstable situations 
where Maxwell-Boltzmann is the limiting form are covered by the extended integrals. 
Extended integral in the standard non-resonant case 

The extended integral in the standard non-resonant case is (Haubold and Kumar, 
2007) 

poo 

ha^ / x^-'[l + a{a-l)x]-^^e-''''~'dx. (4.1) 
Jo 

Theorem 4.1 (Haubold and Kumar, 2007) 

poo 

ha = / x^"'[l + a(Q;- l)x]-^e-''^"'dx 
Jo 

p[a(a - i)]7r (^) ^ ^ i(T,i),(o,i) ; ^ ^ 

where a > 0, 6 > 0, p > 0, a > 1, 3fJ(s) > 0, 3fJ(7 s) > 0. 

When a — )■ 1, ha becomes h- But ha contains all neighborhood solutions for various 

values of a for a > 1. 

If in the above result - is an integer then by taking - = m we obtain 
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Corollary 4.1 . For a > 0, b > and a > 1, we have 

/ x'-^[l + a{a-l)x\-^^e-^'' '"da; 
Jo 



By setting 7 — l = z/, a = l and p = |, we obtain 
Corollary 4.2 . For 6 > 0, a > 1 



/ x''[l + {a-l)x\-—^e-^'''^dx 
Jo 



(o; - i)-+ir (^) '-'v 4 

Extended cut-off case 

In the case of non-resonant thermonuclear reactions with high energy cut-off the ex- 
tended integral is (Haubold and Kumar, 2007) 



fd 

Aa = / x^-\l-a{l-a)x]^i 
Jo 



e-'^-'dx 



where d — ^ , a < 1, a > 0, S — 1, rj — 1, 1 — a(l — a)x > 0, p > 0, & > 0. 



i(l-a) 

Theorem 4.2 



[ x"'-^[l - a(l - Q;)x]T^e-''^ 'dx (4.5) 

" p[a(l-«)]7 ""^^^ 1(7,1), (0,i) ; -^2a i4.bj 

where a > 0, b > 0, p > 0, a < 1, ^{j + s) > and d < 00 

When q; — > 1, 7^^^ becomes 12'^. But contains all neighborhood solutions for various 
values of a for a < 1. In the above result if ^ is an integer then by taking ^ = m we 
obtain 

Corollary 4.3 . For a>0, b > 0, a < 1, d < 00 and 3?(7 s) > 
x'^-^[l - a{l - a)x]^e-'^~^dx 



f 

Jo 



V^(27r)^r + 1) p f a{i - 



[a(l-a)]^ 10, 
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By setting ^ — 1 — a — 1 and p = |, we obtain 

Corollary 4.4 . For 6>0, a<l, d < oo and 3?(i/ + 1 + s) > 

rd 



Jo 



T^i,3 7 lo,i.+i j (4.8) 



Extended depleted case 

Proceeding similarly as in the case of Iia we get for the depleted case 



Theorem 4.3 

poo 

-^""''dx 

1 „9 I / 1 + ? ,1 



/"OO 

ho. = / a;^-^[l + a(a - l)x'']"^e- 



w/iere a > 0, 6 > 0, p > 0, 5 > 0, a > 1, 3fJ(s) > 0, 3f?(7 + s) > 0. 

For the non-resonant thermonuclear reactions with depleted tail ^ — 1 — a — 
1, p = I then we get, 

Corollary 4.5 

/"°° 1-1 
/ + (a- l)a;^]-^e-''^ 



^ r(^) V 



(0,2) 



(4.10) 



By setting f = s' and 25 = m, m = 1, 2, • • • we get 
Corollary 4.6 

2(27r)V^m-i / (a - l)6"V-^+2(l±^ 



ha 



I / 1 N 



( lo ,r,A,..,i^,2(i±;i) ) (4-11) 



4.1 Series representations 

In the following we derive series representations of the right-hand side of (4.4) which 
will be helpful for the evaluation of the extended reaction probability integrals in the 
Maxwell-Boltzmannian form. Taking i/ as a general parameter one can consider several 
situations. Then following through the process in the papers of Haubold and Mathai, 
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see for example Mathai and Haubold (1988), we have the following series representa- 
tions for the extended integral in (4.4). 



4.1.1 Case (I): i/ ±|, A = 0, 1, 2, • • • 

Here we apply residue calculus on the G-function for obtaining the series represen- 
tation of the integrals. Consider the G-function in (4.4). 

-j^ nc+ioo 

1 . \ /(a -1)62^ 



3,1 ^ (a-l)6^ |2-^+ 
^1-3 I 4 lo,i,.+i 



r(.)r(- + . 



xr(i + i/ + s)r 



— V — 1 — s 



ds 



(4.12) 



a-l ' ' J \ 4 

The right hand side is the sum of the residues of the integrand. The poles of the gammas 
in the integral representation in (4.12) follows. 

Poles of r(s) : s = 0, -1, -2, • • • ; T (| + s) : s = -|, -|, • • • ; r(l + i/ + s) : s = 
—V — 1, —V — 2, —V — 3, ■ ■ ■ . 

These arc all simple poles under the conditions in case(l). Then the G-function has a 
simple scries expansion. 
We know that 



lim (s -|- r)r(s) = 



(1 - a)r 
T{a){a)m 

when r(a) is defined, r = 0, 1, 2, • • • ; V (|) = tt^. 



r(a-r) = 
Tia + m) = 



f-l)Tfa) 



(4.13) 

(4.14) 
(4.15) 



a{a + l)---{a + r -I) if r > 1, a 7^ 



1 if r = 0, 

The sum of the residues corresponding to the poles s — —r, r — 0,1,2, ■ ■ ■ is given by 



Ri 



E 

r=0 



r! V2 7 Va-l 



— 1 — u + r ] r(l + u — r) 



TT^r 



^ -i-i/) r(i-i-i/)iF2 ^ 



a — 1 



a — 1 



— 1 — u\ - , —V 
' 2' 



(a - V)\? 



(a - 1)6^ 



(4.16) 

where pFq{ap; hq, z) denotes the generalized hypergeometric function defined as above. 
The sum of the residues corresponding to the poles s = —§,•■• is 



R2 



-27r5r 



a-l 2 \2 



(q - 



X1F2 



1 1 3 1 (a -1)6^ 

-, p; 

a - 1 2 ' 2' 2 ' 4 



(4.17) 



Finally the sum of the residues corresponding to s = —v — 1, —v — 2, —v — 3, • • • is 



^ / 1 3 (a - 1)6^ 



(4.18) 



Adding -Ri, -R2, -R3 we obtain the final result: 

Theorem 4.4 //i/ ^ ±|, A = 0, 1,2, • • • is an integer, then for h > a > 1, we have 



roc 

/ x''[1 + (q;- l)x]"^ 
Jo 



le" 



dx 



r(^-i-'^)r(i + ^) 

^ / 1 1 (a - IW 



2r(^-i--)r(i + .) 



[(« - i)]''+^r (^) 



(a - 1)62 



X1F2 



1 1 3 1 (a -1)^2 

^^~2~''' 2'2~''' 4 



+ 



[(a-l)]^+i 



(a - ly 



XiFsl -;2 + i/,- + z/; 

a — 1 z 



4 

(a - 1)62 



(4.19) 



It is to be noted that the series on the right-hand side of the equation (4.19) are term- 
wise integrable over any finite range. 



4.1.2 Case (II): u is a. positive integer 

In this case some poles of r(s) and r(l + u + s) will coincide with each other. 
Therefore these will be of order 2. We note that the poles s = 0, —1, —2, • • • , —v are 
each of order 1: s = — |, — |, — |, ■ ■ ■ are each of order 1: s = — z/ — 1, — z/ — 2, — z/ — 3, ■ ■ ■ 
are each of order 2. Taking the sum of residues at the poles s = 0, —1, —2, • • • , —u; at 
s — —I, — I, — I, • • • ; at s = —1/ — 1, —u — 2, —1/ — 3, • • • we have 



R.^n'^rf-^^-l-u) r(i + z.)^ 

^ ^ r=0 



(a - 1)6' 



21 



(4.20) 



i?2 



-27r5r 



a- 1 2 



X1F2 



a 



1 
2 



3 1 
' 2'2 



(q - 1)62 



{a - 



(4.21) 
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( 1 \i,2\ 1+^^ °^ 



r-=0 



[a - 1)6' 



2 \ ^ 



In 



+ A 



Br, (4.22) 



where 



and 



^-^ - - + V' + + -0(1 + r) + -0(2 + + r) 

ir^^^^r H - .) r (^) 



(I + z/)^r!(l + z/ + r)! 



From the above results we have the following theorem. 
Theorem 4.5 If i/ > is an integer, then for b > 0, a > 1, we have 



poo 

/ x''[l + {a-l)x]~^ 
Jo 



(tt) 2 



TT^r 



Te-^^"Mx 
1 



[(«-i)r^r(^) 

^(-1)^ (^-1-^) 



X 



a — 1 
a -1)62 



1-1/ r(i + i/) 



-27r2r 
X1F2 

+ 



1 1 



a-1 2 
1 1 



1 



3 1 



I/; -, u- 

a-1 2 ' 2' 2 



(a - 1)6^ 



{a - l)b' 



1)62 y 1)52 xr- 



r=0 



In 



[a - l)b^ 



(4.23) 



(4.24) 



Br } (4.25) 



where Ar and Br are as given in (4-^3) and (4-^4) 



4.1.3 Case (III): ly a negative integer 

Let u — —jjL, n > 0. Then the poles of the G-function in (4.12) are 5 = — —§,••• 
are each of order 1; s = — z/ — 1, —u — 2, — z/ — 3, ■ ■ ■ , +1 are each of order 1; poles 
s = 0,— 1,— 2, ■•• are each of order 2. Again, following through the same stages as 
above we have 
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Theorem 4.6 For u a negative integer, 



(« - ^r^'^ (^) 
X1F2 ^ ^ 



3 1 



a — 1 



2-"' 2'2-"' 



a-l 2 

(a - 1)6^ 



i/ r - + 1/ 



(a - 1)6^ 



+ r(-i-z/)r(---z/)r 



a — 1 



(a - 1)6^ 



-v-2 



X 



E 



■ir (^). 



- r! (1+^)^(2 + ,.), 



(a - 1)6^ 



+ E 



r-=0 



[a - 1)6 



2\ ^ 



where 



— r 



In 



1 



(a - 



5', 



— 1 — + r ) + ■(/'(I + r) + '(/'(r — z/) 



and 



a — 1 

(-i)^+-r(|-r)r(^-i-^-r) 

r!(r — i/ — 1)! 



4.1.4 Case (IV): u a positive half integer 

Let 1/ = m + |, m = 0, 1, 2, • • • . Then 



r(s)r 
r(.)r 



r(i + z/ + s)r 



1 \ 3 

- + s]T{m + - + s)T 



a — 1 

1 

a — 1 



— v — 1 — s 
3 

— m — - — s 



(4.26) 



(4.27) 



(4.28) 



Here the poles of r(| + s) and T{s + m + |) coincides with each other. These will 
be of order 2. We note that the poles s = 0, —1, —2, • • • are each of order 1; s = 
— |, — I — 1, — I — 2, ■ ■ ■ , — I — m are each of order l;s = — m — |, — m — |, • • • are each 
of order 2 and we have the following theorem 
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Theorem 4.7 For v a positive half integer, namely v — m -\- \, m = 0, 1,2, 

, _i 
x''[l + {a-l)x\-—^e-^'' 





■i . , ^ 1 7'"^r — - — m \ T \ - +m 

^a-l)i+-r(^)l \a-l 2 ) \2 

{a - 1)62 



X1F2 



1 1 



a - 1 2 ' 2' 2 



2-^r(^-i-m)r(i + 



(a - l)fe2 



X 



+ 



E 

r=0 



r=0 



Dr\, (4.29) 



where 



Cr = ip{ -m - ^ - r j + ^ ' + ^^(1 + r) + V'(2 + m + r) (4.30) 



and 



(-iy^"""-r(-;»-|)r(^)(^)^ 

(I + m)^r!(l + m + r)! 



(4.31) 



4.1.5 Case (V): v a negative half integer 

Let V — —m — 5, m = 0, 1, 2, • • • . Then in this case 



r(s)r(^ + 5') r{i + u + s)r (-^ 



u — 1 — s 



= rMr(i + .)r(.+ l-™)r(-l^-l + m-. 

Then the poles of Meijer's G-function in (4.12) are s = — r, r = 0, 1, 2, ■ ■ ■ of order 1 
each; s = m — | — r, r = 0, 1, 2, • • • , m — 1 of order 1 each; s = — | — r, r = 0, 1, 2, ■ ■ • of 
order 2 each. Then we have the following theorem 
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Theorem 4.8 For v a negative half integer, namely v — —m — |, m = 0, 1, 2, 



a;''[l + (a-l)x]-^e-''^ 



1 1 



X1F2 



1 1 



m r 



1 



a — i z z z 



a- 1 2 

{a - 1)62 



+ m 



+r ( ?7T. — - ) r(m)r 



a — 1 



X 



E 



r=0 



^' (i-^).(l-^)r 



+ 



r=0 



;« - 1)62 



\a - 1)62 



(a - 1)62 



In 



,2\ 



[a - 1)62 



I^'rk (4.32) 



where 



and 



D' 



+ m + r +V'(l + r)+V'(l + m + r) (4.33) 



(4.34) 



a — 1 

^( -1)"'^^T(^ + ;» + ;■) (^) 
(|)^r!(m + r)! 



4.2 Series representation for the extended cut-off case 

4.2.1 Case(I): 7^ ±|, A = 0, 1, 2, • • • 

Here we apply the same techniques that we apphed previously. Consider the evalu- 
ation of the G- function in (4.8). 



pd ^ 
/ x^ll - {I - a)x\^- 
Jo 



r (t^ + 1) 1 r+'~ r(s)r {s + \)t{i + v + s) 



7r(l - aY+^ 2m 



r(2 + ^+T^ + ^) 



(1 - a)62 



ds(4.35) 



The poles of the integral are as follows: r(s) : s = 0, —1, —2, • • • ; F (| + s) : s — 
-f' -§'■■■ ' r(l + z/ + s) : s = -z/ - 1, -z/ - 2, -z/ - 3, ■ ■ ■ . 
These are all simple poles under case (1). Then evaluating the sum of residues we have 
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Theorem 4.9 //z/y^±f, A = 0,1, 2, • • • is an integer, then for b > 0, a < 1, we have 



1 1 



(1 - a)b^ 



x"[l- (l-a)a;]T^e-*^""da; 

+ / (7r)^r(l + ^) f . 

2(7r)^r(| + i/) f{l-a)b^Y ^ f 1 1 3 1 (1 - q;)^^ 



2' 4 



+ 



r(i + T^) 



^ , 1 3 (1-aW 
1F2 1 + ; 2 + u,- + u; - ^ 



1 — a 



(4.36) 



Here the series on the right-hand side of the equation (4.36) are in computable forms. 

4.2.2 Case (II): u is a. positive integer 

The poles of the gammas in the integral representation of (4.35) are as follows: 
r(s) : s = 0,-l,-2,--- ; r(i + s) : s = -|, -|, • • • ; r(l + + s) : s = 
—u — 1, —v — 2, —u — 3, ■ ■ ■ . 

Some poles of T{s) and r(l + z/ + s) will coincide with each other. Therefore these 
will be of order 2. We note that the poles s = 0, —1, —2, • • • , —u are each of order 1; 
s — — |, — |, — |, • • • are each of order 1; s — —v — 1, —v — 2, —v — 3, • • • are each of 
order 2. Evaluating the sum of the residues we have 
Theorem 4.10 If v > Q is an integer, then for 6 > 0, a <1, we have 



/ a;''[l-(l-a)a;]T^e 
Jo 



, 1 

— -bx'2 



dx 



r(i^ + i) [ (7r)^r(i+^^H-^-i:^) 



a/ r 



v/^(l-a)'^+nr(2 + . + ^)^ (ll(--).H 



(1 - a)b'' 



2(7r)ir(i + z.) ( {l-aW y • 

r(i + ^+I^) V 4 ) ' '\~2 a-1' 2'2 

(l-a)62y+'^g^(l-a)62\'- 

r=0 



1 3 1 il-aW 
u: 



4 



(4.37) 



where 



Er^ip ( -^-ly-rj -'^(^3^ + 1-'^) + -0(1 + +"0(2 + 1^ + (4.38) 



anc? 



;_i)r-(_i)w.r(-i-z.)r(-^) 

(i + ^).H(l + i. + r)!r(^ + l) 



(4.39) 
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4.2.3 Case (III): v a negative integer 

Let u — —/I, II > 0. Then proceeding as in the above cases we have the result 
Theorem 4.11 

rd 



[ x^'il - (1 - a)x]^e-^'' 
Jo 

r(i^ + i) r-2(7r)ir(i + z.) f{l-a)b'Y 



V^(i-«)-+4r(| + z/ + ^) V 4 

^ , 1 13 1 (l-a)62 

X1F2 u ; -.--u: — 

'2 a-12 2 ' 4 



+ 



where 



r=0 



.=0 r\ {1 + 1^)^(2 + 1.), 



(1 - 



in(ii^l+i.; 



(4.40) 



anc? 



= V ( ^ - r ) - ^ (1^ + 2 + i/ - r ) + -0(1 + r) + V'(r - i^) 



(-ir+-r(i-r) 

r!(r-z/-l)!r(j^ + 2 + z/-r) 



(4.41) 



(4.42) 



4.2.4 Case (IV): u a positive half integer 

Proceeding as before, we have 

Theorem 4.12 For 1/ a positive half integer, namely v — m + \, m = 0, 1, 2, 



x^[l — (1 — (y.)x\ i-«e" 



-63;" 3 



dx 



r(i^ + i) / (7r)ir(m + |) 



-m 



^/^(l-a)'^+H^(m+| + ^) 

2(7r)ir(l + m) / (l-a)6^ y y. (-l)^'(-l 
r (2 + m + ^) V 4 ) ^ r-\ ( 



3 11 

2 ~ 1 -a' 2' 



m 



1 (1 - q;)62 



m 



2' 4 



(l-a)62y"+5^/(l-a)6^ 

r=0 



r=0 
2 \ ^ 



(4.43) 
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where 



Gr = ^l!{ -m -^-A + ^ U + ~ ^ 1 + + r) + ^(2 + m + r) (4.44) 



and 



-m 



I) (---i), 



r!(l + m + r)!r(l + ^) (l + ^) 



(4.45) 



4.2.5 Case (V): v a negative half integer 

Let u — —m — |, m = 0, 1, 2, • • • .In this case we have 

Theorem 4.13 For v a negative half integer, namely v — —m — |, m — 0,1,2, 



x^ll - (1 - a)x] i-^'e 



dx 



r(i^ + i) f (vr)ir(i + m) / 1 



1 1 1 

2 ^ 1 - a' 2' 2 



m; 



(1 - a)62 



V^(l-a)-+ilr(|-m+^) 

^ r(m-i)r(l + m) |^ (i-a)6^ y-gi (-i^), f{i-a)b 



r(i + T^) 



2\ 



r=0 



(1 - a)b' 



2\ 5+'- 



where 



and 



H' 



(-ir(-^-T^). 

r!(m + r)!r (l + — m) 



(4.46) 



= ^ ( - r j + ^ f 1 + ^j-^^ m-r] + ^(1 + r) + ^(1 + + r) (4.47) 



(4.48) 



5 Behaviour of the integrals ha and /. 



id) 
2a 



The behaviour of the integral Iia is such that as the value of the pathway parameter 
a changes the curve will move away from the stable situation ie, the strict Maxwell- 
Boltzmann situation(Figure 3 below). The graphs of the integral Iia when u — 1 and 
at a = 1, a = 1.25, a = 1.35, a = 1.45 are plotted in Figure 1. As a > 1.5 the 
G-function nolonger exists as it violates the conditions. We can take other value of v 
also. 
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"(d) 



v=l 

(a) a=l 

(b) a=1.25 

(c) a=1.35 

(d) a=1.45 




10 12 14 IB 18 2Q 



Figure 1. Behaviour of Iia for various values of a > 1 

Similarly the behaviour of the integrals Jg^'' is such that the function moves away from 
the stable case and comes closer to the origin. 



.)a=l 
(b)a=0.75 
c) a=0.65 
(d) 0=0.55 




Figure 2. Behaviour of for various values of a < 1 

As a — )■ 1 we get the standard situation which is done in the series of papers of Mathai 
and Haubold. As a — ?■ 1 the two integrals will come close to the following limiting 
situation. 



v=l 5=1 

(a) a=1.35 — 

(b) a=1.25 

(c) a=l — 

(d) a=0.75 ■■■■ 

(e) a=0.65 — 




0.4 



D.3 



D.2 



0.1 



a=l 
v=l 



10 12 14 IB 



10 20 30 40 50 

(a) (b) 
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Figure 3. Maxwell-Boltzmann case or the limiting situation a = 1. 

In figure 4 as the value of 6 moves we can see the depletion in the high energy tail of 
the Maxwell-Boltzmann situation. The graphs of depletion in the stable situation as 
well as many unstable and chaotic situations are plotted here. The cases (i),(ii) and 
(iii) show the depletion when a = 1.25, a = 1.35 and a = 1.45 respectively and (iv) 
shows the depletion in the stable situation (a = 1). 




(iii) (iv) 
Figure 4. Depletion for 5 = 1, 2, 3 and a = 1.25, 1.35, 1.45, 1 

6 Conclusion 

For the analytic evaluation of the probability integral for equilibrium conditions we 
consider a more general form of the reaction probability integral. We investigated in 
section 4 the series representations of the extended integral, Iia and Jg^'' where as the 
pathway parameter a — )■ 1 one gets the Maxwell-Boltzmann case. Iia for various values 
of a is plotted in Figure 1, for various values of a is plotted in Figure 2 and in Figure 
3 the limiting case,namely a = 1 case or Maxwell-Boltzmann situation is plotted. The 
plotting is done by Maple 9. 
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